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Abstract 

The quantum phase transitions of metals have been extensively studied in the rare-earth "heavy 
electron" materials, the cuprates, and related compounds 1 . The Fermi surface of the metal often 
has different shapes in the states well away from the critical point. It has been proposed^* 3 - that 
these differences can persist up to the critical point, setting up a discontinuous Fermi surface change 
across a continuous quantum transition. We study square lattice antiferromagnets undergoing a 
continuous transition from a Neel state to a valence bond solid, and examine the fate of a small 
density of holes in the two phases. Fermi surfaces of charge e, spin 1/2 quasiparticles appear in 
both phases, enclosing the usual Luttinger area. However, additional quantum numbers cause the 
area enclosed by each Fermi surface pocket to jump by a factor of 2 across the transition in the limit 
of small density. We demonstrate that the electronic spectrum across this transition is described 
by a critical theory of a localized impurity coupled to a 2+1 dimensional conformal field theory. 
This critical theory also controls the more complex Fermi surface crossover at fixed density, which 
likely involves intermediate phases with exotic fractionalized excitations. We suggest that such 
theories control the electronic spectrum in the pseudogap phase of the cuprates 4 . 
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I. INTRODUCTION 



Many strong correlation problems of current interest relate to the physics of doped Mott 
insulators. Most observed metallic states can be understood in conventional terms: ordering 
associated with broken spin or lattice symmetries, along with Fermi surfaces of charge e, 
spin 1/2 quasiparticles. Nevertheless, the details of the state differ depending upon the 
importance we attach to strong correlation physics e.g. different Fermi surface topologies 
and magnetic moment configurations are obtained from the weak and strong interaction 
limits, namely (i) a spin density wave instability of the free electron band structure, or (u) 
a spin-wave theory of a lightly doped magnetically ordered Mott insulator. It is clear that 
rich physics lies in understanding the interpolation between such limiting regimes, and that 
there a number of important applications e.g. to the widely discussed differences between 
the underdoped and overdoped cuprates. 

A central question is whether there could be a single continuous quantum phase tran- 
sition between such states with distinct types of order or Fermi surface topologies (recent 
experiments^ on CeRhIn 5 are compatible with an abrupt or very rapid change in Fermi sur- 
face topology). Such critical points could be the key to understanding the many observed 
'non-Fermi liquid' properties. However, a complete critical field theory of a continuous tran- 
sition between such states has so far been presented only in insulators^, where the two states 
have broken symmetries which cannot be connected by a continuous transition in Landau 
theory; instead a 'deconfined' theory focussing not on order parameters but on fractionalized 
excitations and emergent gauge forces shows that such transitions are indeed possible. An 
extension of these ideas to metallic systems has been discussed^, including a scenario which 
would allow for a discontinuous change in a Fermi surface. Others^ have proposed a 'local 
quantum critical' scenario for transitions between states with different Fermi surface topolo- 
gies, based upon numerical studies in extensions to dynamic mean field theory. However, 
the analytic structure of this mean field transition remains unclear, and it is not clear how, 
even in principle, this 'local critical' structure can be extended beyond mean field theory. 

In this paper, we will build on the theory of the deconfined quantum critical point pro- 
posed in Ref. [6 for insulating spin S = 1/2 antiferromagnets on the square lattice^^ 1 ^ 1 ^; 
the transition is between (see Fig. [1]) a state with magnetic Neel order (with spin polarized in 
the staggered configuration) to a spin-gap state with valence bond solid (VBS) order (which 
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FIG. 1: Schematic phase diagram. The ellipses represent spin singlet valence bonds. The coupling s 
tunes the insulator across the Neel- VBS transition, and S is the mobile hole density. The deconfined 
quantum critical point is at s = s c in the insulator with 5 = 0. The vacancies ('holons') carry a 
gauge charge q = ±1 under an emergent U(l) gauge force. In the cartoons above, the reader can 
interpret q as a sublattice label. In the s < s c , Neel phase, q determines the spin: a vacancy on an 
up (down) spin site carries net spin down (up) and so is equivalent to a charge e spin- 1/2 hole. For 
s > s c , the hole is a composite of a vacancy and a nearby unpaired spin with opposite q, moving 
by rearranging nearest-neighbor valence bonds; note that this motion preserves spin and sublattice 
quantum numbers separately (see also Ref. |9|). So there are twice as many states per momentum 
for a charge e spin 1/2 hole in the VBS state than there are in the Neel state. 
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FIG. 2: Momentum space Fermi surfaces in the Neel and VBS regions of Fig[TJ The filled circles are 
the 4 K p wavevectors, with K\ = (ir/2a)(l, 1), K 2 = {ir/2a){\, -1), K s = -K\, K4 = -K 2 , with 
a the lattice spacing. The dashed line in the Neel phase indicates the boundary of the magnetic 
Brillouin zone. Only the Fermi surfaces within this zone contribute to the Luttinger counting, and 
so the area of each ellipse is Af = (27r) 2 <5/4. In the VBS phase, all 4 pockets are inequivalent, 
and so the area of each ellipse is Af = (27r) 2 <5/8. The dashed lines now show the reduction of 
the Brillouin zone due to the VBS order which appears at sufficiently low temperatures; "shadow" 
Fermi sufaces, with weak photoemission intensity (estimated in the text), will appear as reflections 
across these lines, and these Fermi surfaces are not shown. 
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is spin rotation invariant, but breaks lattice symmetries by ordering of valence bonds). We 
will describe the metallic states obtained by doping this system with a small density (5) of 
holes. Figs. [TJ and [2] summarize our key results, and the captions support them with simple 
physical arguments. The decoupling of the spin and U(l) gauge charge in the VBS phase 
described in Fig. [TJ and the consequent doubling of the hole flavors in Fig. [21 leads to our 
precise result: the discontinuity in lim^o Af / 'S on the two sides of the Neel-VBS quantum 
critical point, as shown in Fig. [21 

On the other hand, the evolution of the Fermi surface along the fixed 5 line P in Fig [1] 
involves separate questions associated with the interplay between the mobile holes and the 
s = s c deconfined critical point. If the 6 = monopole- induced confinement physics^ survives 
for S > 0, then there will be a single critical line which is crossed by P at which the Fermi 
surface changes discontinuously. However, this appears unlikely because the monopoles can 
be screened by the charge carriers^, leading to a narrow intermediate fractionalized 'holon 
metal' phase with no conventional Fermi surface within the shaded region. In either case, 
there is a wide momentum, frequency, and temperature regime over which the hole spectral 
function is incoherent. We will show here that this is described by a universal theory of 
a single localized vacancy^ coupled to the 2+1 dimensional conformal field theory of the 
5 = 0, s = s c critical point. The incoherent, critical hole spectral function is momentum- 
independent, and so while this view of the transition has a 'local' character-, the degrees 
of freedom include a non-trivial, bulk quantum field theory. Further, we believe that the 
evolution of the fermionic excitations across the shaded region of Fig [1] cannot be captured 
by dynamic mean field theory methods or self-consistent impurity models. 



II. FIELD THEORY 



The motion of charge carriers doped into an insulating square-lattice quantum antiferro- 
magnet is conventionally described by the "t- J" model, 

H t -j = - Uj(clc ja + h.c.) + Y J iA • Sj, (1) 

where c ia creates an electron with spin a on the sites i of a square lattice and Si = 
I Ylaa' c i(T&<r<7' c iv'- m addition, the constraint c\ a Ci a < 1 is enforced on each site, mod- 
elling the large local repulsion between the electrons. It is important to note that our results 
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are more general than a particular t-J model, and follow almost completely from symmetry 
considerations. 

We begin our analysis by recalling the physics in the absence of doped holes; we then 
only need to consider the second term in Eq. ([T]) i. e. the "Heisenberg model" . With weakly 
frustrated Jjj the ground state is the Neel state. With stronger frustration (or additional 
ring exchanges), there can be a quantum transition to a VBS phased. The field theory for 
the Neel- VBS transition in the insulator 6 is expressed in terms of a complex 'relativistic' 
field z a , where a =|, j, which is the field operator for a charge neutral, spin S = 1/2 spinon 
excitation. Low lying singlet excitations are represented by an emergent U(l) gauge field 
A^, where // = r, x, y is a spacetime index. The field theory has the action J d 2 rdrC z with 

A = \(d» - iA^)z a \ 2 + s\z a \ 2 + | (\z a \ 2 ) 2 + ... (2) 

Here s is the tuning parameter in Fig. [TJ we have set a spinon velocity to unity, and u is a 
spinon self-interaction. For s > s c , this field theory is in a spin SU(2)-invariant U(l) spin 
liquid phase, with gapped spinon excitations, and provides a continuum description of the 
Schwinger boson mean field theory of the square lattice antiferromagnet^. The inclusion of 
monopole tunneling events shows that this spin liquid is ultimately unstable to confinement 
and the development of VBS order—. This monopole physics is peripheral to our main 
results, and so we defer its discussion until later. 

Now we need to add charge carriers to this antiferro magnet. We have represented the spin 
excitations of the antiferromagnet by Schwinger bosons above, and so the charge carriers 
require a spinless fermion /, the 'holon' (see Appendix [A] for more details). Such a mean- 
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We can take the 



field lattice study has been carried out some time ago in Refs. 
low energy limit of their results, and so obtain the field theory needed for our analysis. 
However, we can also proceed on symmetry grounds, when our main ingredients are the 
transformation properties of the various fields under square lattice symmetry operations. 
Because the elementary quanta also carry U(l) gauge charges, it is not required that the 
fields be invariant under such operations, only their gauge-invariant combinations. This 
larger space of transformations was dubbed the projective symmetry group (PSG) by Wen2£, 
and we list the PSG of the z a in Tabled! This / holon also carries the gauge charge of the U(l) 
gauge field in Eq. (j2J). Indeed, because the spinon field flips charges under a sublattice 
interchange in Table [J we see that holons on the two sublattices have charges q = ±1. This q 
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TABLE I: PSG transformations of the fields under square lattice symmetry operations. T x : transla- 
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tion by one lattice spacing along the x direction; R d ^: 90° rotation about a dual lattice site on the 



plaquette center [x — > y,y —> —x); Z£ ual : reflection about the dual lattice y axis (x — > — x,y — ► y); 
T: time-reversal, defined in real time. e a/ g is the antisymmetric tensor with eti = 1 and e a p = e a ^ . 
The local VBS order parameter is defined as ^vbsIt) = ( — l) r ' x {S? ■ S? +x ) + £(— l) r ' y (S? ■ SV+y); 
where x,y are the unit lattice vectors. 

quantum number will play a central role. However, the holons can also acquire an additional 
'valley' quantum number if their dispersion minimum is not at the center of the Brillouin 
zone; we will keep track of this at the cost of some additional bookkeeping. The analysis of 
the holon Green's function of the present U(l) spin liquid phase, along the lines of Ref. Il9 . 
shows that there are 2 distinct valleys, with labels v — 1, 2, and we choose a gauge in which 
their minima are at wavevectors K i 2 shown in Fig. [2j Summarizing, we have 4 species of 
spinless fermions f qv , with U(l) charge q = ±1, and valleys v = 1,2. All the information 
we need to construct the effective action and the observable correlations of these fermionic 
holons is encapsulated in their PSG, which is derived in Appendix |XJ and is summarized in 
TableE 

An immediate consequence of the PSG of the z a and the f qv is that we can write down 
an expression for the physical charge — e, spin-1/2 electron annihilation operator c a (f) = 
X/p=i e lKp ' r ^ 'paif) where the K p are defined in Fig. [2] and 



#1,3* 
#2,4a 



z af+l i £a/3 z ^ f—l 
z af+2 ^ £ a/3 z ^ f-2 



(3) 
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These relations will allow us to construct the Fermi surfaces in Fig. [2j 

The field theory of the doped antiferromagnet has the Lagrangian is £ = C z + £/ + C c . 
The fermionic holon terms are 

= E & - ^ " ^ - ( ^2m^ j)2 ) ^' (4) 

where j extends over x, y, m m = and m^x = vrt\y are the mass of the elliptical hole 
pockets, fi is the hole chemical potential, and the x and y directions are rotated by 45° from 
the principle square axes. The PSG prohibits a linear derivative term in Eq. (j4j), and so 
the f qv dispersion is an extremum at zero momentum; this pins the Fermi surfaces to be 
centered at K p in the Neel phase (Fig. [2]), but not, as we will see, in the VBS phase. Finally, 
we include coupling between the holons and spinons invariant under all PSG operations. 
There is an unimportant scalar coupling ^ \z a \ 2 fq V fqv, but the more important term is 

C c = i\e a(3 (flif-iZadsZp + fl 2 f-2Z a dyZ^ + cc, (5) 

equivalent to the dipole coupling introduced by Shraiman and Siggia^, arising from the 
hopping of electrons between nearest-neighbor sites. Another symmetry-based derivation 
of the theory of the doped antiferromagnet was recently giverJ 2 . for the Neel phase alone. 
In contrast, our analysis starting from a spin rotation invariant spin liquid can address 
the quantum critical point, and can also be specialized to the Neel state by spontaneously 
breaking spin rotation symmetry. We now turn to an analysis of the motion of holes in the 
field theory £ = C z + Cf + C c , first in the Neel and VBS phases, and then finally at the 
deconfined quantum critical point that seperates these two phases. 

III. DOPING THE NEEL STATE 

First, consider the s < s c Neel phase of C z with (z a ) 7^ 0. By the Higgs mechanism, 
the condensate of z a renders the A M massive, and so we can safely integrate the out. A 
crucial point is that the Higgs mechanism also ties the U(l) gauge charge q of the holons 
to the spin quantum number S z along the Neel order (because of the broken symmetry, S x 
and S y are not conserved); more precisely S z = q/2, and so the f qv carry all the quantum 
numbers of the electron. This is established in Appendix 02 
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Now we consider the Fermi liquid state obtained with a total density S of the four fermion 
species f qv . Each will form a separate Fermi surface containing 5/4 states; their coupling 
to the spin-wave can be treated perturbatively, and do not modify any of the key Fermi 
liquid characteristics. From Eq. ([3]) we deduce that there are 4 hole pockets centered at the 
K p wavevectors, each enclosing the area Af = (2vr) 2 (5/4), as shown in Fig. [2J The caption 
of Fig. [2] shows that the same area is obtained in direct Fermi liquid counting of electrons 
within the magnetic Brillouin zone. 



IV. DOPING THE VBS STATE 



Next, we examine doping the s > s c state. For now, we neglect monopoles, and so the 
undoped state is a spin rotation invariant U(l) spin liquid with representing a gapless 
'photon' excitation; the spinful excitations are gapped, but single S = 1/2 spinons do not 
exist as asymptotic states — the logarithmic 'electrostatic' potential mediated by the A T 
binds the spinons in pairs. We can identify 2A as the gap towards creating the lowest 
spinful excitation; this vanishes at the quantum critical point as A ~ \s — s c \ u . In a theory 
with iV spinon species (at — 1 ... N in Eq. (T5])), the logarithmic potential will have the form 
V(r) = (12A/JV) ln(rA) in the large N limit^. Now add a single holon into this state. 
The coupling to A T causes this hole to have an 'electrostatic' self energy which diverges 
logarithmically with system size^, and so a spin-singlet single holon state is not stable. 
Rather, the holon will peel off a single spinon from above spinon gap, and form a S = 1/2, 
charge e bound stated. This bound state has the same quantum numbers of as an ordinary 
hole, and is neutral under the A^ U(l) gauge force. A finite density of such bound states 
can then form a Fermi surface with charge e S = 1/2 quasiparticles. We also have to 
consider the pairing between holons with opposite U(l) gauge charges, induced by the A^ 
gauge force (this attractive force must be balanced against the repulsive Coulomb force 
associated with the physical electromagnetic charge e of each holon). This is expected to 
lead to superconductivity at low enough temperatures, but we will not discuss this here; our 
focus is on the normal state. 

The structure of the holon-spinon bound state is described in more detail in Appendix O 
As a consequence of the mixing between the different bound states induced by the terms 
in C c , the bound state eigenmodes turn out to be precisely the electron-like momentum 
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eigenstates in Eq. with a dispersion which has a minimum at K% = ±A where A ~ AA 2 . 
This is responsible for the shift in the centers of the elliptical hole pockets shown in Fig. [21 

The usual counting argument now allows us to deduce the area of each hole pocket. There 
are 4 inequivalent pockets, and a factor of 2 degeneracy for spin; so Af = (27r) 2 (5/8). The 
factor of 2 difference from the result obtained in the Neel phase is one of our key results. 

The state we have described so far is actually not a conventional Fermi liquid: the area 
enclosed by its Fermi surface is not the same as that of a non-interacting electron gas at 
the same filling and with the same size of unit cell. In the non-interacting case one would 
have a Fermi surface area (27r) 2 (l — 5)/2. In our case, where the holes were doped on a 
background spin liquid state with a gapless photon excitation, we obtain that the electron- 
like area enclosed by the Fermi surface is (27r) 2 (l — 5/2) . This state is a fractionalized 
Fermi liquids, obtained here in a single band model, in contrast to its previous appearance 
in Kondo lattice models. 

The instability of the spin liquid to confinement and VBS order induced by monopoles^ 1 ^, 
and the associated halving of the Brillouin zone will finally transform the state into one 
obeying the conventional Luttinger theorem, with an electron-like area enclosed by the 
Fermi surface of (27r) 2 (l — S)/2. The VBS order parameter is a complex field, $vbs ; which 
is proportional to the monopole creation operator. Its expectation value vanishes upon 
approaching the critical point as ($vbs) ~ A |0VBS / 1 '. The exponent /?vbs/^ is the scaling 
dimension of the monopole operator, which is expected to be large^ 2 ^. We are interested 
here in the mixing induced by ($vbs) across the reduced Brillouin zone shown in Fig [21 
However, because of the shift in the minimum of the holon-spinon bound state dispersion 
discussed above, this mixing is negligible as long as y/5 < A ~ A 2 , and so can always be 
neglected as 5 — > 0. This establishes the Fermi surface structure claimed in Fig. [2J 

At larger 5, the hole pockets will cross the dashed lines in Fig. [2, and Bragg reflections 
will split the Fermi surfaces. We analyze this magnitude of the mixing matrix element in 
the Hamiltonian, using the PSG of $vbs m Table 1, in Appendix [D] We find that it has 
a value ~ A($vbs) ~ A 1+/3vBS ^. This splitting of the Fermi surface is negligible provided 
the matrix element is smaller than the hole kinetic energy, or 5 > A 1+/3vBS ^. Over such 
a possible regime of larger 5, the basic pictures of Figs. [H and [2] continue to hold. The 
photoemission intensity of the "shadow" Fermi surfaces noted in Fig. [2] is proportional to 
the square of the matrix element or ~ A 2 ( 1+/3vBS /^ 
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V. QUANTUM CRITIC ALITY AND CONCLUSIONS 



To complete the picture, let us address the physics of the shaded region in Fig. HJ Here 
we have a theory of a finite density of holons f qv interacting with spinons z a via a U(l) gauge 
force Ap. At sufficiently long scales, the holons will 'Thomas- Fermi' screen that longitudinal 
A T force, and so obviate the binding into gauge neutral combination. Consequently the 
spinons and holons remain as relatively well-defined excitations, and we enter a fractionalized 
holon metal phase. This screening can be prevented if the spacing between the holons 
(~ is larger than the holon-spin binding length ~ 1/A, or 5 < A 2 ; this fixes the 

boundary of the shaded region in Fig. [TJ The shaded region will exclude the unshaded 
region of the VBS phase where the Fermi surface splitting is negligible (discussed in the 
previous paragraph) provided /?vbs > u , an inequality that holds at least for large N. 

We now describe the criticality of the hole spectrum at s = s c . We assume here an 
observation scale (frequency (u) or temperature (T)) large enough, or a 5 is small enough, 
so that the holes can be considered one at a time. A key observation about a single hole is 
that its quadratic dispersion (i.e. the terms proportional to l/(2m V x,y) in Cf) is an irrelevant 
perturbation on the quantum critical point of C z which involves excitations which disperse 
linearly with momentum. Consequently, the hole may be considered localized, and itsphysics 
is closely related to the single impurity in a spin liquid problem analyzed in Ref. [l5|. Here 
we are interested in the single hole Green's function and this requires the overlap between 
states of the spin liquid with and without the impurity. This can be computed by analyzing 
the quantum critical theory of C z coupled to one holon localized at r = 0, represented by 
the r-independent Grassmanian fij): 

S qc = J d 2 rdrC z + J drf (J- + s - iA T (r = 0, T )j f. 

Here Eq is an arbitrary energy fixing the bottom of the holon band, and, following earlier 
arguments^, the only relevant coupling between the 'impurity' holon degree of freedom and 
the bulk degrees of freedom of C z is the gauge coupling associated with A T . The spectral 
function of a physical charge e, S = 1/2 hole is given by the two-point correlation of the 
composite operator z a f* . If this operator has scaling dimension rjh/2, then the universal 
critical hole Green's function, Gh, is independent of wavevector and obeys the scaling form 

G h = T-^^(h(cu-s )/k B T) (6) 
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where $ is a universal scaling function. At T = 0, we obtain an incoherent spectrum 
associated with the power-law singularity Gh ~ (uj — eo)~ 1+Vh . We have computed rjh by a 
standard 1/N expansion for Gh under the action S qc and found 



Another perspective on the exponent rjh is obtained by mapping it to an observable in 
the lattice non-compact CP 1 (NCCP) model studied by Motrunich and Vishwanath 13 . On 
a three-dimensional cubic lattice with spacetime points j, the model has the complex spinor 
fields Zj a on each lattice site, and a vector potential A, M on each link extending along the 
\i direction from site j. The hole Green's function at imaginary time M T (in units of the 
lattice spacing) is then given by the two-point correlator of the Zj a in the time direction 
along with an intermediate 'Wilson line' operator (representing the contribution of the / 
fermion) which renders the correlator gauge invariant (see also Ref . [2^ ) : 



this is expected to decay as e~ E ° T /MJ/ h for large M T at the quantum critical point. The 
continuum limit of the above correlator was computed by Kleinert and Schakel^ in the 1/N 
expansion: their result agrees with Eq. (JTj). 

This paper has presented a theory for a strongly interacting quantum critical point in 
a metal. The strong interactions imply scaling of observables as a function of thj/ksT, 
which is often seen experimentally, but is absent in the weak-coupling theory^. The theory 
has a interesting evolution in the Fermi surface geometry across the transition, which is 
reminiscent of that found in the solution of self-consistent impurity models^. However, the 
details are different, and can only be addressed in a quantum field-theoretic framework which 
we have described. This work therefore also serves to place studies of quantum transitions 
in the 'dynamical mean field theory' approach^>21 in the field-theoretic context. 

The T dependent broadening in the quantum-critical hole spectrum in Eq. ([6]) has sim- 
ilarities to recent observations in the underdoped cuprates 4 . The dispersion minimum is 
shifted away from the K p points in Fig. [2} as is the case experimentally, and in contrast to 
other theories of the normal state based upon Dirac fermion spinons. This suggests that a 
theory such as ours, based upon a fractionalized Fermi liquid stated, should lead to a useful 
description of the pseudogap phase. 




(7) 
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APPENDIX A: DERIVATION OF THE PSG 

We are interested here in the physics of the t-J like models in which the Hilbert space is 
restricted to 3 possible states on each site of the square lattice: an up/down spin electron, 
and a vacancy. The state with two electrons is projected out because of the strong on-site 
repulsion, and this constitutes the strong correlation physics. We can represent these states 
in terms of canonical operators by decomposing the electron operator into spinon and holon 
operators. On one sublattice of the square lattice we write the electron operator, c a as 

c a = b a fl (Al) 

where b a are canonical 'Schwinger' bosons and /+ are canonical fermionic holons. The 
projection onto the 3 allowed states is imposed by the local constraint 

flh + b a K = l (A2) 

on each site. On the other sublattice, we use bosons which transform as a conjugate respre- 
sentation 

c a = e a pb f! (A3) 

with a similar constraint. 

Now consider the mean-field U(l) spin liquid state of the insulator. This is described by 
the effective Hamiltonian^^ 

hj = -qJ2 b ^ + H - c - + A E 6 * 6 * + A E *S<S ( A4 ) 

(ij) i 3 

where i is restricted to be on one sublattice, and j on the other, and Q, X are positive 
constants. Any operation which interchanges the two sublattices leaves Hj invariant under 
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the PSG mapping 



b a — > e a pb 

T -> e a %. (A5) 



Further, the invariances of Eqs. flAl|A3l) under such an operation demand that 

/+ - f- 

f- - -/+■ (A6) 



Note especially the sign in the last equation, which is a consequence of £ Q/3 £« 7 = —5". 

The spinon spectrum of Hj shows that the minimum energy excitations are near zero 
momentum, and so we may take its low energy limit by a naive gradient expansion. This 
leads^ to the effective theory C z in Eq. (j2J) with 

z a ~b a + P a (A7) 

The PSG of the z a in Table U follows from the above relations. 

Determination of the PSG of the f qv requires the additional information that the holons 
have their dispersion minima K\i. The phase factors associated with the transformation 
properties of e lKv ' r , combined with those in Eq. (1A6[) lead to the results in Table fl] 



APPENDIX B: HOLES IN THE NEEL STATE 



By spin rotation invariance, we can always rotate the z a condensate (and without any 
rotation of the spinless f qv ) to produce a Neel order N a = z a *<j^fzp (where a a are the 
Pauli matrices) polarized along the (0, 0, 1) direction. Now examine the response of the 
theory to a uniform magnetic field H applied along the z direction. Under such a field, 
the only change in the action is that 15 (d T — iA T )z a — > (d T — iA T + (H/2)a z )z a . Choosing 
(z a ) = y/\s — s c \/ud a ^ we obtain a term in the Lagrangian 

£ = • • • (\s - s c \/u) [iA„ - (H/2)5 flT ] 2 + ■ ■ ■ (Bl) 

which gaps the photon. Integrating out the A^ and then evaluating M z = 5S/5H\ H=Q , 
we obtain an expression for the magnetization density M z 

M z =±Y,<lflU + --- (B2) 

qv 
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where the ellipses represent an additional term which measures the magnetization of the 
spin waves. This establishes, as claimed in Fig. [Tj, that S z = q/2 for the fermions in the Neel 
phase. 



APPENDIX C: HOLON-SPINON BOUND STATES 

We describe the structure of a holon-spinon bound state in more detail. The bound state 
of a f + i holon with a spinon will be mixed by C c with the bound state of a /_i holon with 
a spinon (parallel considerations apply to the f±2 holons). 

Before writing the wave-function for this bound state, we need to decompose the rel- 
ativists field z a into non-relativistic, canonical boson fields p a (the spinon) and h a (the 
anti-spinon). At low momenta, this is done by the standard parameterization 

z a = -JL= ( Pa + e a(3 h m ) (CI) 

In the absence of C c the bound state of f + i and h a and the bound state of /_i and p a 
are independent, and thus the problem can be solved in the familiar center-of-mass frame, 

M- = m vl + A; Rj = + ^ (C2) 

P] = -r~\] r = r h -r s (C3) 

A is the mass of the spinon in the non-relativistic limit, and m v -j are the holon masses in Cf. 
The holon co-ordinate is represented by fh and that of the spinon by f s . The kets describing 
the motion of the bound state can be labeled in terms of the momentum P of the bound 
state center-of-mass, 



\ap+) = j dr-sdrle^^fUf^ifsM 

\aP-) = [ dr- s dr^ P M^fUr h )p a Hr s )\0) (C4) 



p2 

The kets correspond to the eigenvalues = ^fe and have been written in terms of the "first- 
quantized" ground state wavefunction <f>(r) of the holon-spinon interaction, V(r). V(r) = 
(12A/JV) ln(rA) binds the holons and spinons over a length scale ~ 1/A, with 0(0) ~ A. 

The mixing of the ± kets due to L c can be studied by interpreting this term simply as 
a Hamiltonian 7i c with the operator identification in Eq. (IC1I) . 7i c only connects kets with 
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the same P and a with the matrix element (aP + \H c \aP—) = — P * (we used the fact 
that d x (j)(0) = 0). Now, by simply diagonalizing a 2 x 2 matrix we can infer the dispersion 
Ep induced by C c to first order in perturbation theory in A, 



h ? ~ 2Mj ± M x [L5} 

with eigenmodes which correspond precisely to the electron states in Eq. ([3]) (x is replaced 
by y for f±z) . So the bound state dispersion has a minimum at K% = ±A where A = 
A 1 0(0) | 2 ~ A 2 . This is responsible for the shift in the centers of the elliptical hole pockets 
shown in Fig 2. 

APPENDIX D: COUPLING TO VBS ORDER 

Here we describe the terms in action which couple the VBS order parameter $vbs to 
the holons and/or spinons. Such terms are responsible for the mixing between the Fermi 
surfaces of the VBS state (shown in Fig. [2]) and the 'shadow' Fermi surfaces obtained by 
Bragg reflection across the reduced Brillouin zone boundaries (shown as dashed lines). 

The most general terms in the action follow by the requirements that they be U(l) gauge 
invariant, and also invariant under the PSG in Table 1. These requirements turn out to be 
extremely restrictive. 

First, we search for terms which involve $vbs an d a bilinear of the f qv fermions. However 

a detailed analysis shows that there no such terms which are invariant under all the PSG 

operations; this is seen by first listing the f qv bilinear invariants under /^ ual (under which 

<3>vbs is invariant), and then noting that their transformations under R^f are incompatible 

with those of $vbs- Consequently, the coupling between the VBS order and the fermions 

will be weaker than might have been initially expected, and will vanish faster than ($vbs) ~ 
A/Svbs/* as A __, o. 

The simplest non-vanishing coupling turns out to require the full ty p electron operators 
in Eq. ([3]). This has the form 



£vBS — A VBS ^ VBS 



- ^l^i + \&^ 3 - ^3^2) 



+ (^ 2 - + #J*3 - ^4)1 + CC. (Dl) 
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Now we need an estimate of the matrix element of £vbs between the holon-spinon bound 
states found in the previous section. Using Eqs. (I CI I) and flC4j) , we obtain a matrix element 
of order 

Avbs($vbs>^^ ~ A VB sA 1+/3vB ^. (D2) 

As expected, this does vanish faster than A /3vBS / iy , and is responsible for the weak Bragg 
reflection across the reduced Brillouin zone boundaries of the VBS state in Fig. [2j 



1 P. Coleman and A. J. Schofield, Nature 433, 226 (2005). 

2 Q. Si, |cond-mat/0601001| 

3 T. Senthil, S. Sachdev, and M. Vojta, Physica B 359-361, 9 (2005); T. Senthil, Annals of 
Physics, 321, 1669 (2006). 

4 A. Kanigel, M. R. Norman, M. Randeria, U. Chatterjee, S. Suoma, A. Kaminski, H. M. Fretwell, 
S. Rosenkranz, M. Shi, T. Sato, T. Takahashi, Z. Z. Li, H. Raffy, K. Kadowaki, D. Hinks, 
L. Ozyuzer, and J. C. Campuzano, Nature Physics 2, 447 (2006). 

5 H. Shishido, R. Settai, H. Harima, and Y. Onuki, J. Phys. Soc. Jpn. 74, 1103 (2005). 

6 T. Senthil, A. Vishwanath, L. Balents, S. Sachdev, and M. P. A. Fisher, Science 303, 1490 
(2004); T. Senthil, L. Balents, S. Sachdev, A. Vishwanath, and M. P. A. Fisher, Phys. Rev. B 
70, 144407 (2004). 

7 T. Senthil, S. Sachdev, and M. Vojta, Phys. Rev. Lett. 90, 216403 (2003); T. Senthil, M. Vojta, 
and S. Sachdev, Phys. Rev. B 69, 035111 (2004). 

8 Q. Si, S. Rabello, K. Ingersent, and J. L. Smith, Nature 413, 804 (2001). 

9 M. Levin and T. Senthil, Phys. Rev. B 70, 220403 (2004). 

10 A. Kuklov, N. Prokof'ev, B. Svistunov, and M. Troyer, Annals of Physics 321, 1602 (2006) have 
studied a particular lattice realization of the 'easy-plane' limit of the field theory proposed in 
Ref. (j, and found a first-order transition. This does not rule out a second-order transition for 
other lattice realizations. A second-order, deconfined, transition has been theoertically estab- 
lished for N sufficiently large, but finite. There is also strong numerical evidence for a continuous 
transition for the physically interesting case with N = 2 and full SU(2) symmetr y 11 ' 12 ' 13 . 



11 A. W. Sandvik, [cond-mat/0611343 

12 M. Kamal and G. Murthy, Phys. Rev. Lett. 71, 1911 (1993). 



17 



13 O. I. Motrunich and A. Vishwanath, Phys. Rev. B 70, 075104 (2004). 

14 M. Hermele, T. Senthil, M. P. A. Fisher, P. A. Lee, N. Nagaosa, and X.-G. Wen, Phys. Rev. B 
70, 214437 (2004). 

15 A. Kolezhuk, S. Sachdev, R. R. Biswas, and P. Chen, Phys. Rev. B 74, 165114 (2006). 

16 D. P. Arovas and A. Auerbach, Phys. Rev. B 38, 316 (1988). 

17 N. Read and S. Sachdev, Phys. Rev. Lett. 62, 1694 (1989). 

18 C. L. Kane, P. A. Lee, and N. Read, Phys. Rev. B 39, 6880 (1989). 

19 Y. M. Li, D. N. Sheng, Z. B. Su, and L. Yu, Phys. Rev. B 45, 5428 (1992). 

20 X.-G. Wen, Phys. Rev. B 65, 165113 (2002). 

21 B. I. Shraiman and E. D. Siggia, Phys. Rev. Lett. 61, 468 (1988). 

22 F. Kampfer, M. Moser, and U.-J. Wiese, Nucl. Phys. B 729, 317 (2005); C. Briigger, F. Kampfer, 
M. Moser, M. Pepe, and U.-J. Wiese, |cond-mat/0606766[ 

23 G. Murthy and S. Sachdev, Nucl. Phys. B 344, 557 (1990). 



24 W. Rantner and X.-G. Wen, |cond-mat/0210152] L. 

25 H. Kleinert and A. M. J. Schakel, Phys. Rev. Lett. 90, 097001 (2003). 

26 M. T. Beal-Monod and K. Maki, Phys. Rev. Lett. 34, 1461 (1975); J. A. Hertz, Phys. Rev. B 
14, 1165 (1976); A. J. Millis, Phys. Rev. B, 48, 7183 (1993). 

27 A. Georges, G. Kotliar, W. Krauth, and M. Rozenberg, Rev. Mod. Phys. 68, 13 (1996). 



18 



